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Abstract. - In a two-dimensional microwave chaotic cavity ohmic losses located at the 
contour of the cavity result in different broadenings of different modes. We provide an analytic 
description and establish the link between such an inhomogeneous damping and the complex 
(non-real) character of biorthogonal wave functions. This substantiates the corresponding 
recent experimental findings of Barthelemy el al. [Europhys. Lett. 70 (2005) f62]. 



Open wave-chaotic systems in the presence of energy losses (absorption) are nowadays 
under intense experimental and theoretical investigations, see [1,2] for recent reviews as well 
as [3] for a general discussion. Most of the works concern the case of uniform absorption which 
is responsible for homogeneous broadening Thorn of all the modes (resonance states). However, 
in some experimentally relevant situations like, e.g., complex reverberant structures [4,5] or 
even microwave cavities at room temperature [6,7] one should take into account also localized- 
in-space losses which lead to an inhomogeneous part Tj n h of the widths which varies from 
mode to mode. As a result, the neighboring modes experience nontrivial correlations due 
to interference via one and the same decaying / dissipative environment that result in the 
complex-valued wave functions of corresponding resonance states. Such a complexness may 
reveal itself in long-range correlations of wave function intensity and current density [8] that 
were recently studied experimentally [9]. Following Refs. [4, 10], it is convenient to measure 
the above mentioned complexness through a single statistical parameter, namely the ratio 
((lmi[)) 2 ) I ((Reip) 2 ) — q 2 of variances of the real and imaginary parts of the mode wave function 
ip. The modes are real (q = 0) in the case of vanishing inhomogeneous losses and become 
complex-valued when T; n h ^ (the value of Thorn has no effect on the mode complexity, as 
it will become clear later on) . The strong experimental evidences in the favor of the intimate 
relation between q and Ti n h were recently provided by Barthelemy et al. [6], who analyzed 
hundreds of resonances in a Sinai-like chaotic microwave cavity and found a proportionality 
between these two quantities. Here, we present a thorough analytic description of this problem. 

© EDP Sciences 
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It is instructive to start with the discussion of the physical picture of the problematic. 
A typical experimental setup consists of a flat (two-dimensional) cavity fed with microwaves 
through attached antennas or waveguides which support M a scattering channels (propagating 
modes) in total. These very channels are used to perform all the measurements and their 
number is finite. Dissipation through ohmic losses at cavity boundaries gives another reason 
to treat our system as open, since boundaries may be viewed locally as distributed "parasitic" 
equivalent channels with ad hoc impedances [11,12]. One should distinguish, however, be- 
tween almost uniform ohmic attenuation at the cavity plates and localized absorption at the 
cavity contour [6,7]. The number of "parasitic" channels responsible for the former ("bulk") 
mechanism can be naturally estimated as M b ~ (L/A) 2 , while one has typically M c ~ L/X 
channels at the contour (where L is the characteristic length of the cavity and A denotes the 
wave length). Both Mb,M c » 1 but their ratio M c /Mb ~ X/L is parametrically small and 
that will be essential for our consideration. 

It is natural, therefore, to use the following model description of the problem. According to 
the Hamiltonian approach to scattering [13-15], see also [2], one can represent the scattering 
matrix in terms of the effective non-Hermitian Hamiltonian H. e s of the open system as follows: 

s tot = i-iyt ^_ 1 v> n eS = h - l -vvK (i) 

The Hamiltonian H of the closed chaotic system gives rise to N levels (eigenfrequencies) e„ 
characterized locally in the relevant range of the energy E by the mean level spacing A. 
Those are coupled to all the open channels via the N x (M a + Mb + M c ) matrix V of coupling 
amplitudes and, as a result, are converted into complex resonances £ n = E n — |T„, which are 
given by the poles of the S'-matrix. Accordingly, we propose to decompose V — {A, B, C} 
into coupling to M a antennas, Mb "bulk" and M c "contour" parasitic channels to separate 
explicitly different contributions to the widths T n . The total (M a + Mb + M c )-dimensional 
scattering matrix (1) is, of course, unitary( 1 ). However, one can access experimentally only 
the (M a x M a ) subblock 5=1- iA^(E - H c g)~ 1 A which is subunitary. 

Without loss of generality one can consider the case of preserved time-reversal invariance 
(which was indeed the case of experiment [6]) when the coupling amplitudes V are real and 
H is symmetric. It is convenient first to represent 7i e ff in the eigenbasis of its Hcrmitian part 
as follows: 

(M a M b m c \ 

X! A n A m + ^2 B n B m + ^2 C n C rn j ■ ( 2 ) 
a=l b=l c=l / 

In a chaotic cavity, H is commonly described in the framework of Random Matrix Theory [3] . 
A (real orthogonal) rotation that diagonalizes the random Hcrmitian matrix H transforms the 
(fixed) matrix V to Gaussian-distributed coupling amplitudes with zero means and following 
covariances (we assume statistical independence of channels from different classes) [13,14]: 

(«',) - 2 Ka -5 aa 'S nn , , (B b n B«) = 2n b -5 bb '8 nnl , (C£C#> = 2k c -S cc ' '<$„»< • (3) 

7T 7T 7T 

Coupling constants k determine transmission coefficients T = 4k/(1 + k) 2 of the corresponding 
channels, so that T< 1 (T = 1) stands for weak (perfect) coupling. The strong inequality 
Mb 3> M c 3> 1 allows us to perform now the limit of a very large number of weak fictitious bulk 
channels, M b — » oo and T b — > with M b T b = 27rr hom /A being kept fixed [12], which singles 

( 1 )Unitarity of the total S-matrix as well as causality, T n > 0, are automatically provided by the factorized 
algebraic structure of the anti-Hermitian part of Heff [14]. 
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out the homogeneous absorption contribution. Indeed, by virtue of the central limit theorem 
one may replace in 7i e s the sum J2b=i ^n^m w i tn lis average value (^n^m) = Thom^nm 

in the limit considered that yields 

(H cS ) nm = [e n - l -T hom ^j 5 nm - l - (AA T + CC T ) nm , (4) 

meaning that all the levels acquire one and the same attenuation rate. Since 7i c ff comes into 
the scattering problem only as the resolvent (E — Hcs)^ 1 , uniform absorption turns out to be 
operationally equivalent to a pure imaginary shift of the scattering energy E — > E + |rh om = 
Ej [12, 16]. Thus the physical scattering matrix acquires the following form: 

S = l- iA T 1 -^A , H cH = H- 1 -{AA T + CC T ) . (5) 

E 1 — H c g 2 

A representation similar to (5) was used in [5, 17] to study statistics of transmitted power in 
dissipative ergodic microstructures. Here, we concentrate rather on spectroscopic problems. 

It is clear from the above consideration that only escape to antennas and inhomogeneous 
losses contribute to the fluctuating part T n = T n — Thom of the ^widths. The latter are 
given now by the imaginary parts of the complex eigenvalues of H e ff. Since H c e is non- 
Hermitian, the eigenvalue problem H c fi\n) = £ n \n) and (n\H c s = (n\£ n defines two sets of 
right and left eigenfunctions, which satisfy the conditions of biorthogonality, (n\m) — S nm , 

and completeness, X^=i l n )("-l = Ijv- The matrix U nm = (n\m) ^ S nm differs from the unit 
one and is known in nuclear physics as Bell-Steinbcrger nonorthogonality matrix [18] (see a 
compact description in [14]). This matrix features in two-point correlations in open systems 
seen, e.g., in decay laws [19]. U nn appears also in optics via the so-called Petermann factor of 
a lasing mode [20-22] . Some statistical aspects of chaotic nonorthogonal eigenfunctions were 
recently studied in [22-24]. 

For a general non-Hermitian matrix left and right eigenvectors are independent of each 
others. However, in our case 7i e ff is complex symmetric (due to time-reversal invariance), 
implying that left eigenfunctions are given by the transpose of the right ones, (n| = |n) T . As 
a result, TC c s = ^£ ^ T can be diagonalized by a complex orthogonal transformation [14], with 
£ = diag(£i, • • • , £jy) and ^ = (|1), • • • , \N)), that leads to the well-known pole representation 
of the iS-matrix, 5=1 — iA T &(E 1 — £)~ ll & T A 7 or in its components: 

Saa' - S aa , -if^^L , C = A a \n) =f^Al\n) k . (6) 

n=l ^1 k=1 

The wave function component of the n-th mode excited through the a-th channel is gen- 
erally complex. This complexness is solely due to biorthogonal nature of the eigenfunctions 
and is directly related to the structure of the anti-Hermitian part of H e s- In particular, all 
■0^ would be real were (luiH c e)nm oc S nm or, more generally, if the anti-Hermitian part of the 
effective Hamiltonian commuted with its Hermitian part( 2 ). 

We proceed now with considering the case of tunneling coupling to antennas, which was 
realized in experiment [6]. This allows us to neglect safely antenna contributions to T„, 
approximating H c s ~ H — \CC T . In the case of the large but finite number M c of (contour) 



( 2 )In this case both the Hermitian H and anti-Hermitian parts of H e g can be diagonalized simultaneously, 
thus the eigenbasis being a conventional orthogonal one, as follows from H T = H. 
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channels, the levels acquire on average the width given by the so-called Weisskopf's estimate 

r inh = M C T C A , (7) 

well-known in nuclear physics, see, e.g., [14,15]. It is worth noting that this value can be 
formally linked to Sabine's law of room acoustics, which determines the average width r re fl 
related to absorption at the cavity contour. One has r re fl = cLT c /(irS) [6,25], where S is the 
cavity area, c is the speed of light and L is now the cavity perimeter. Making use of Weyl's 
law for the mean level spacing A = cX/S and putting M c = L/(A/2), we find that r re fl is 
exactly converted to T m h. This provides us with a further link between the present model 
description and the microscopic treatment of [6,7] based on Maxwell's equations. 

Fluctuations of the widths around Ti n h are mostly due to those of the matrix CC T . As 
follows from the central limit theorem, fluctuations of off-diagonal matrix elements are sup- 
pressed as compared to diagonal ones at M c » 1, (CC T )„^ m <~ (CC T ) nn /^M c <~ T- m h/ \/M c , 
so that they contribute to T n in the next-to-leading order in 1/M C . However, off-diagonal 
matrix elements give a dominating contribution to the mode nonorthogonality. Indeed, all 
essential features of the problem can be most explicitly seen in the two-state approximation. 
Representing H e s as follows 

HcS ~ { e 2 ) 2{ (Cjd) ||C 2 || 2 ) ' (8) 

where are the M c -dimensional vectors of the corresponding coupling amplitudes, one can 
easily solve the secular equation for its two eigenvalues £1,2, finding E\.2 = (ei+e2 z fd)/2, where 
d = \J (li - e 2 ) 2 - (C T C 2 ) 2 and 61,2 = £1,2 - \ 1 1 Ci , 2 1 1 2 - The corresponding eigenfunctions (we 
assume ei < £2) are given by 

m-«(! f ) ^ /=^L, t.) 



Af 2 = (1 — ,/ 2 ) -1 being the normalization constant. Then the nonorthogonality matrix reads 
U = |W| 2 [(1 + |/| 2 )l2 + 2Re(f)a y ], with the Pauli matrix a y . The parameter / controls the 
mode complexness, as follows from 

ff i2 = #(^ 2 ipy (10) 



for the mode wave functions. 

In experiment [6], the complexness parameter 

q 2 = ((ImC) 2 )/((ReC) 2 > (11) 

was accessible only in the regime of the weak mode overlap due to inhomogeneous losses (we 
stress, however, that the total width Thom + rinh can be larger than A). In this regime, r; n h <C 
A, the computation of q can be easily carried out by making the use of the above two-state 
approximation. One finds straightforwardly that q 2 w ((Re/) 2 ) w ((C T C 2 ) 2 )/(2A) 2 , which 
determines now the average nonorhogonality matrix as (U nm ) w (1 + 2q 2 )S nm . The remaining 
Gaussian averaging over random amplitudes (3) yields J2 C c'(CiC%Cf C% ) = M c (2k c A/it) 2 w 
M c (T c A/2tt) 2 , as T c w 4k c <C 1 in our case of large M c and fixed Ti n h. Collecting everything, 
we arrive at 

1 Tinh „„s 
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that substantiates the proportionality between q and Ti n h established experimentally [6] . The 
mode complexness (nonorthogonality) decreases as the number of absorptive channels in- 
creases, in agreement with the general discussion presented above. 

We discuss now the role of fluctuations in energy levels and contributions from the other 
levels neglected so far. Restricting ourselves to the same weak overlap regime, we can use 
the perturbation theory for wave functions to find q 2 w j^2 m ^ n ((CjC m ) 2 )((e n — e m )~ 2 ). 
The known correlation function i?2(w) of Gaussian orthogonal ensembles, see [3], can be now 
used to get X^m#n(( e « — £ m)~ 2 ) = J duJLd~ 2 R.2(u>) , which can be appreciated as (average 
square of) the so-called "level curvature" studied in [26,27]. Since i?2(w) ~ w as w — > 0, 
this integral has a logarithmic divergency regularized by setting the lower limit ~ that 
also ensures us to stay within the perturbation theory. As a result, (12) is renormalized to 
yield a contribution ~ ln(A/r; n h). One should expect that such a factor would be absent in 
a non-perturbative treatment (see also the relevant discussion in [22]). The computation of q 
at arbitrary inhomogeneous absorption is still an interesting open problem to consider. 

In summary, we have presented the model description for inhomogeneous (localized-in- 
space) losses in open chaotic systems and discussed thoroughly the resulting complexness 
(or biorthogonality) of the mode wave functions. In particular, the complexness parameter q 
determining the relative weight of the imaginary parts of the modes is found analytically to be 
proportional to the inhomogeneous part r; n h of the widths in full agreement with experimental 
results of [6]. Though the present calculation of q is perturbative in the small parameter 
r inh /A, it may nevertheless be valid in the intermediate or large modal overlap regime where 
Thom dominates, as, for instance, in the case of room acoustics or elastodynamics [1,25]. Our 
analysis may be also relevant for problems of mode nonorthogonality outside of the scattering 
systems as those considered recently in [23,28]. 
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